Evolution of the Quantum Friedmann Universe Featuring 

Radiation g 



O 



O 

(N 



V. V. Kuzmichev 

Bogolyubov Institute for Theoretical Physics, National Academy of Sciences of Ukraine, 

Metrolohichna St. l^b, Kiev, 03143 Ukraine 

Abstract: The classical and quantum models of the Friedmann universe originally 
filled with a scalar field and radiation have been studied. The radiation has been used 



to specify a reference frame that makes it possible to remove ambiguities in choosing the 



time coordinate. Solutions to the Einstein and Schrodinger equations have been studied 
under the assumption that the rate of scalar-field variation is much less than the rate of 
universe expansion (contraction). It has been shown that, under certain conditions, the 



quantum universe can be in quasistationary states. The probability that the universe 

goes over to states with large quantum numbers owing to the interaction of the scalar 

and gravitational fields is nonzero. It has been shown that, in the lowest state, the scale 

factor is on order of the Planck length. The matter- and radiation-energy densities in 

^ ■ the Planck era have been computed. The possible scenarios of Universe evolution are 

O ■ discussed. 
(N 

§ ' 1. INTRODUCTION 

O ' 

That quantum gravity theory cannot rely on experimental data [1] adds importance to 
exactly soluble cosmological models. However, the application of basic ideas underlying 



£jrj! quantum theory to a system of gravitational and matter fields runs into difficulties of a 

fundamental character, which do not depend on the choice of a specific model. By way 
of example, we will consider a homogeneous, isotropic, and closed universe characterized 
by the Friedmann- Robert son- Walker metric; that is, 



als 2 = a 2 (ri)[N 2 (ri)dr] 2 -dn 2 ], (1) 

Here, N(rj) is a function that specifies the time-reference scale; a(r]) is a scale factor; dfl 2 
is an interval element on a unit 3-sphere; and rj is the parameter that is related to the 
synchronous proper time t by the differential equation dt = Nadrj. Considering that 
scalar fields play a fundamental role both in quantum field theory and in the cosmology 
of the early Universe [2, 3], we assume that, originally, the Universe was filled with matter 
in the form of a uniform scalar field <p. If the field <p varies slowly in the early Universe, 
its potential V(<f>) specifies the vacuum-energy density (cosmological term) and ensures 
Hubble expansion. Restricting our analysis to the case of minimal coupling between 
geometry and the scalar field, we represent the action functional in the conventional 
form 

S= [ d V [^a' + Tt^'-H], (2) 
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where a prime denotes differentiation with respect to d/drj; 7c a and tt^ are the momenta 
canonically conjugate with the variables a and 0, respectively; and H is the Hamiltonian 
given by 



H =-N 
2 



a 



= NTZ. (3) 

Here, the variables a and are taken, respectively, in units of the length I = a/2 G/3tt 
and in units of = a/3/8 ttG. The function N plays the role of a Lagrange multiplier, 
and the variation 5S/5N leads to the constraint equation TZ — 0. The structure of 
the constraint is such that true dynamical degrees of freedom cannot be singed out 
explicitly. This creates problems in the interpretation of quantum geometrodynamics 
[46]. It is commonly thought that the main reason behind such difficulties is that there 
is no natural way to define a spacetime event in general covariant theories [7]. In the 
model being considered, the above difficulties are reflected in that the choice of the time 
variable is ambiguous. 

For the choice of the time coordinate to be unambiguous, the model must be supple- 
mented with a coordinate condition. When the coordinate condition is added to the field 
equations, their solution can be found for a fixed time variable. However, this method 
of removing ambiguities in specifying the time variable does not solve the problem of 
a quantum description, because undesirable consequences of this ambiguity in eventual 
equations cannot be avoided in this way. 

In this study, we propose specifying a reference frame with the aid of an additional 
matter source. This method does not come into conflict with the adopted ideas of the 
early Universe [2, 3]. At the same time, it enables us to study the evolution of the 
Universe not only in the semiclassical approximation but also at a purely quantum level. 

2. CLASSICAL DESCRIPTION 
2.1. Fundamentals of the Model 

The ambiguity associated with choosing the time coordinate in (1) will be removed 
with the aid of a coordinate condition imposed prior to varying the action functional, but 
its coordinate invariance will be restored [7, 8]. We will choose the coordinate condition 
in the form T' = N, where T is the privileged time coordinate, and include it in the 
action functional with the aid of a Lagrange multiplier P; that is, 

S = Jdr ] [n a a' + 7r <p( j ) ' + PT'-n}, (4) 

where 

n = N[p + n] (5) 

is the new Hamiltonian. The constraint equation reduces to the form 

P + TZ = 0. (6) 

Integrating the canonical equation P' = [P, Ti] = 0, we immediately obtain P = E, where 
E is a constant. The full set of equations for the model in question becomes 

d 2 -y0 2 + t/ = £, (7) 
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<i ■ 2 dV 



+ 2- + a A — 

a d<p 



0, 



(8) 



where overdots denote differentiation with respect toT and U = a? — 4 V(0). Equation (7) 
represents the Einstein equation for the (g) component, while equation (8) is the equation 
of motion for the field 0. A modification to the Einstein equations that is associated with 
including the coordinate condition in the action functional is that, on the right-hand 
side, there additionally arises an energy-momentum tensor T ° = E/a 4 , T 1 = Tf = Tf = 
— E/3a 4 , Tp = for a ^ (5 that can be interpreted as the energy-momentum tensor 
of radiation [9]. The choice of radiation as the matter reference frame is natural for the 
case in which relativistic matter (electromagnetic radiation, neutrino radiation, etc.) is 
dominant at the early stage of Universe evolution. If our Universe were described by the 
model specified by equation (4), it would be possible to relate the above radiation at the 
present era to cosmic microwave background radiation. 

2.2. Solving the Einstein Equations 

A feature peculiar to the model in question is that it involves a barrier in the variable 
a. This barrier, described by the function U, is formed by the interaction of the scalar 
and gravitational fields. It exists for any form of the scalar-field potential V((f>) and 
becomes impenetrable on the side of small a in the limit V — > 0. In just the same way as 
in inflation models (see [2, 3]), we assume that the rate at which the scalar field changes 
is much smaller than the rate of universe evolution, \a/a\ ^> |0|. In this case, equation (7) 
can be integrated in a general form. Presented immediately below are explicit solutions 
in the regions a < a\ and a > a 2 , where they are assumed to satisfy the boundary 
conditions a(0) = and a(ti n ) = a 2 , respectively; here, a\ and a 2 are the turning points 
(di < a 2 ) specified by the condition U = e, and t{ n is the initial instant of time in the 
second region. We have 



a(t) 



^— (l - cosh 2Wt) + J ^ sinh2v / Ut 



1/2 



(9) 



for a < cii and 



a(t) = ^ 217 T 1 + v / l-4Uecosh2v / U(t - t in ) 



1/2 



(10) 



for a > a 2 . The quantities e and U depend parametrically on 0. In the zero-order 
approximation, the former is given by e = E. The above solution to equation (7) can be 
refined by taking into account a slow variation of the field with the aid of the equation 



a 
~2 



+ e = E, 



(11) 



where e(0) stands for a potential term, which is bounded by the inequality E < e < 1/4V. 
The case of e > 1/4 V, which corresponds to an infinite motion, will not be considered 
in this study. From equations (8) and (11), it follows that, in general, a change in the 
potential V(<f>) entails a change in the quantity e(0). 

The solution in (9) describes the universe expanding from the point of the initial 
cosmological singularity to the maximum possible value of ai achieved at the instant 
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t m = —K= In ( after that, the expansion rives way to contraction, and the 




universe collapses by the instant t c — 2 t m . For 2 y/V t <C 1, the solution in (9) takes the 
form 



a(t) ~ [2^t] 



1/2 



(12) 



It is independent of V and describes the evolution of the universe that is dominated by 
radiation [9] and which expands in the de Sitter mode from the point a = a 2 . In the 
extreme case of e = 0, where there is no radiation, the region a < a x contracts to the 
point a = 0, and the expansion can proceed only from the point a = a 2 . Since the region 
a < a 2 cannot be treated in terms of classical theory, it is assumed that the classical 
spacetime with a = a 2 is formed as the result of a tunnel transition from "nothing" 
taken to mean some quantum state of the protouniverse (see, for example, [2, 3,10,11]). 
If, originally, the universe was filled not only with matter but also with radiation, it 
can undergo evolution in the region a < a\ as well. In the general theory of relativity, 
the solutions in (9) and (10) describe two independent scenarios of the evolution. The 
inclusion of the mechanism of quantum tunneling through the barrier U requires a joint 
analysis of these scenarios. It is then legitimate to consider the probabilities of finding 
the universe in each of the classically accessible regions. 

The evolution of the universe depends on the initial distribution of the classical field 
and its subsequent behavior as a function of time. The chaotic- inflation scenario [3], 
which is realized in the region a > a 2 , is described by equations (7) and (8) as applied 
to the case specified by the inequalities (^f^) 2 < 1, V > |^ - ^-|, and jz\(f>\ < |^|. 
In the model where the scalar-field potential V is taken to be proportional to 0™, the 
chaotic-inflation process proceeds between scalar-field values greatly exceeding a level of 
(initial stage) and those achieving this level (final stage). In this approach, radiation 
has virtually no effect on the degree of inflation, and the scalar field represents the field 
of an inflaton. The de Sitter regime of inflation persists as long as the potential V(<f)(t)) 
varies rather slowly with time. From equations (8) and (11), it follows that the inequality 
V + e/a 4 < holds in the expanding universe (a > 0). If the potential V increases with 
time, the quantity e is bound to decrease. But if V decreases, e can increase, and the 
rate of this increase is higher for greater a. We will now estimate e by using the relation 
e ~ TqO 4 . In our Universe, with a ~ 10 28 cm, the main contribution to the radiation- 
energy density comes from cosmic microwave background radiation with energy density 
p° ~ 10" 10 GeV/cm 3 . Setting T ° = p°, we find that, at the present era, the result is 
e = e 7 ~ 10 117 . In the early Universe, the scale parameter is a ~ 10~ 33 cm, while the 
energy density Tq is on the order of the Planck value. On this basis, it can be found that 
e ~ 1 corresponds to that era. It follows that e increased in the evolution process. This 
increase can be explained by a considerable redistribution of energy between the scalar 
field and radiation at the initial stage of Universe existence. Quantum theory is able to 
account for this phenomenon in a natural way (see below). In the region a > a 2 , the 
possible variation of e with time does not affect the quasiexponential expansion of the 
universe, because the inflation stage terminates in a rather short time interval of t ~ 10~ 37 
s [3] , and the evolution process is then determined by other factors (particle production, 
heating, etc.). In the region a < a±, the role of the increase in e with decreasing V may 
prove substantial. In principle, the dependence of e on <p(t) makes it possible to provide 
the missing power in the t dependence of a and to solve the problem of the Universe 
size. In order to demonstrate this explicitly, we assume that, up to the present time 
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to, the Universe has expanded according to the law specified by (12) [2]. We then have 

/ \l/2 

a(to)/a(t p ) = ( A/eo/ep t /t p ) , where eo = e(0(to)), e P = e (0(^p)), an d t p is the Planck 
time. The ratio of e and e p can be estimated as e /e p ~ V p /V , where V p = V((f>(t p )) 
and Vq = V(0(to))- Assuming that, in the Planck era, V p is on the order of the Plank 
energy density and that Vq is on the order of the mean matter-energy density at the 
present era, po = lCT 5 GeV/cm 3 , we find that the value of a(to) ~ 10 28 cm corresponds 
to a(t p ) ~ 10" 33 cm. 

3. QUANTIZATION 
3.1. Schro dinger Equation 

In quantum theory, the constraint equation (6) comes to be a constraint on the wave 
function that describes the universe filled with a scalar field and radiation. Replacing the 
canonically conjugate variables involved in equation (7) by the operators a = ax, 7r a = 
— id a , (j) = 0x, = —id<f,, andP = —idr, we find that the state vector (a,0|\l/(T)) 
satisfies the equation (a, (f)\^(T)} 



2id T \m{T)) 



dl--„dl-U 



a a 2 



|¥(T)). (13) 



where the order parameter is assumed to be zero [5, 10-13]. Equation (13) represents an 
analog of the Schrodinger equation with a Hamiltonian independent of the time variable 
T. The momentum P associated with radiation appears linearly in equation (13). We 
can introduce a positive definite scalar product < oo and specify the norm of a 

state [8, 11]. This makes it possible to define a Hilbert space of physical states and to 
construct quantum mechanics for the universe model being considered. 
A partial solution to equation (13) has the form 

|*(T)> = |^)expj^(T-T )}, (14) 

where the state ip satisfies the time-independent equation 

-d 2 a + ^dl + U-E)j |^) = 0. (15) 

The quantity E is arbitrary in the general theory of relativity, but, in quantum theory, 
it is quantized in accordance with solutions to equation (15). 

3.2. Quasistationary States 

In considering the quantum case, we assume that, at the initial stage, the motions 
occurring in the system under study can be separated into two types: a slow variation 
of the scalar field, in which case the operator (2/a 2 )<9| can be treated as a perturbation, 
and fast changes in geometry. This assumption is a quantum analog of the adiabaticity 
hypothesis, which leads, in the zero-order approximation, to the solutions given by (9) and 
(10). In quantum theory, the problem being considered reduces to solving the equation 

[dt-U + e n (<P)]\vn) = 0. (16) 
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The wave functions ip n and the eigenvalues e n , which depend on parametrically, describe 
the evolution of the universe for very slow variations of the potential V associated with 
the field (more specifically, under the condition \d\nV/d(fi\ <C 1). In order to take into 
account the variations of the field 0, we can represent ip as an expansion in terms of 
the states f n and integrate then equation (15). The quantum number n of the system 
unperturbed by the operator (2/a 2 )d^ will be a good quantum number for the universe 
in the perturbed state ip as well. 

We will further consider solutions to equation (16), allowing for the possible bound- 
ary conditions. For e n < 1/4 V the classically accessible regions a < a x and a > a 2 
are bounded by the turning points a x and a 2 , which are now dependent on the state of 
the quantum system. In the region a > a 2 , a general solution has the form of a super- 
position of converging and diverging waves. In the Wentzel-Kramers-Brillouin (WKB) 
approximation, we can write 

where C\ is the amplitude of an " incident" wave describing the universe whose scale factor 
decreases, while C 2 is the amplitude of the wave "traveling" toward greater values of a 
and describing the expanding universe. For the extreme case of e n = 0, which corresponds 
to the radiation-free universe with an undetermined time variable, the wave function in 
the form (17) was studied by many authors (see, for example, [3, 5, 10-13]). It coincides 
with the Vilenkin wave function [10] at C± — and generalizes the Hartle-Hawking wave 
function [11] to the case of e n ^ 0. 

If we consider a universe formed at an instant separated by a comparatively large time 
interval [— oo < (T — To) < 0] from the commencement of observation, the boundary 
condition C\ — will imply that the diverging wave corresponding to a quasistationary 
state [14] is singled out from the superposition in (17). No situation that is physically 
realizable can exactly correspond to the requirement C\ — for all instants of time be- 
cause, in that case, the process that leads to the formation of a quasistationary state 
and which involves converging waves would be eliminated from the analysis. According 
to the general concepts of quantum theory [14, 15], a quasistationary state can be im- 
plemented approximately by requiring that the region where the asymptotic form (17) 
with C\ = is realized be bounded by the condition a < a max , in which cases y? n is set 
to zero for a > a max , a-max ~ y^T being some boundary value of the scale factor. That 
the instants of time that satisfy the conditions T < T Q and T > a max /^/e^ are excluded 
from the analysis is physically justified because this makes it possible to avoid specula- 
tions about the properties of the scalar field and radiation in the Universe at times that 
are practically inaccessible to observation and for which there are no reliable hints from 
high-energy physics. 

The possibility of quantum tunneling through the region a\ < a < a 2 of the potential 
barrier results in that stationary states cannot be realized in the region a < a±. If, 
however, the potential V(4>) is sufficiently small, quasistationary states whose lifetime is 
much greater than the Planck time can exist in the region a < a±. The probability r„ 
of the decay of the universe occurring in a given quasistationary state ip n can be found 
by requiring that the wave function in (17) satisfy the radiation condition C\ = 0, which 
selects the discrete complex values e„ = e„ + ir„[14]. 
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Table 1: Values of e n as computed in the WKB approximation on the basis of (19) and 
within perturbation theory (PT) for various values of the potential V (here, U max = 1/AV 
is the barrier height; Aa = a2 — a\ Is the barrier width; and X is the number of levels 
at given V), along with the decay-probability values T n as given by (18) under the same 
conditions 



n 


V 


e„(WKB) 


en(PT) 


r„ 


Umax 


Aa 


£ 





0.08 


2.62 


2.63 


0.31 


3.125 


1.03 


1 




0.05 


2.79 


2.79 


0.006 


5 


2.25 


1 




0.03 


2.89 


2.88 


2 xl0~ 6 


8.33 


3.70 


2 




0.02 


2.93 


2.92 


7 xlO- 11 


12.5 


5.08 


3 




0.01 


2.97 


2.96 


1Q -24 


25 


8.05 


6 


1 


0.03 


6.34 


6.35 


0.01 


8.33 


2.06 


2 




0.02 


6.59 


6.59 


10~ 6 


12.5 


3.70 


3 




0.01 


6.80 


6.80 


io- 19 


25 


6.92 


6 


2 


0.02 


9.88 


9.94 


0.003 


12.5 


2.34 


3 




0.01 


10.51 


10.51 


io- 15 


25 


5.93 
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We further impose the boundary condition tp n \ a= o = on the wave function tp n . For 
the case of T n <^ e n , we then obtain 



r„ = 2 



ai I (22 \ 

^ -jJ== expi -2 J y/U-e n da I , (18) 

.0 I ai J 



where e n is determined from the equation 



ai 

J V e n - U da — 7r + . 



(19) 



In the extreme case of V — 0, equation (19) can easily be integrated, which yields 
e n \v=o = ^ = 4n + 3; that is, we can see that, for all values of n, coincides with the 
energy of an isotropic oscillator with zero orbital angular momentum [16, 17]. According 
to (19), the first level (that at e ~ 3) emerges at V ~ 0.08. 

At small V, equation (19) leads to e n values coincident with those that are obtained 
directly from equation (16) by perturbation theory in a 4 V(<f)). The table displays e n 
values calculated by perturbation theory and in the WKB approximation [that is, with 
the aid of equation (19)]. Also presented in this table are the decay-probability values 
r n computed for various potentials V. Since T n Ree n , the decay-probability values 
r n found on the basis of (18) are expected to be close to true values for small n as well. 

We note that the smaller the value of a 4 V at given e n , the higher and the broader 
is the potential barrier U and, hence, the smaller is the decay probability T n . If some 
state ip n is characterized by a small value of r n , the possibility that this state decays 
can be disregarded over the decay time r = 1/T„, so that this state can be considered to 
be stationary in this limit. This corresponds to defining a quasistationary state as that 
which takes the place of a stationary state when the probability of its decay becomes 
nonzero [14]. 
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In describing the universe on the basis of equation (15), the process of universe pro- 
duction from "nothing" in the radiation-free model (E = 0) [2, 10, 12, 13] is replaced 
by quantum tunneling from a quasistationary state with a definite (complex) value of 
E. It can easily be seen that the solution given by (17) describes the de Sitter regime of 
expansion according to (10). In order to demonstrate this explicitly, we note that, in the 
WKB approximation, we have — id a ip n ~ — y/e n — U y? n ; that is, the classical momen- 
tum is given by n a = — a = — y/e n — U, whence we obtain equation (7) in approximation 
\a/a\ > |0|. 

In quantum models not featuring radiation, the De Witt- Wheeler equation [5, 18] 
determines the time-independent wave function of the universe. This leads to well-known 
difficulties in interpreting this function and in comparing results obtained on its basis 
with the observed evolution of our Universe [3, 6]. If, however, the case of E = is 
considered as the limit to which the model with a privileged reference frame reduces 
when E — > 0, we can also speak about the time evolution of the Universe free from 
radiation [19]. 

3.3. Dynamics in the Prebarrier Region 

By studying inflationary scenarios of the evolution of a universe filled with a scalar 
field, it was revealed that a "realistic" potential V must decrease with time [10]. With 
decreasing V, the number of quantum states in which the universe can occur increases, 
while the decay probabilities decrease sharply (see table). The first instants of the ex- 
istence of the universe are especially favorable for its tunneling through the potential 
barrier U. 

A quasistationary state tp n takes the place of the stationary state whose wave function 
in the region a < ai; is close to the wave function \n) of the state unperturbed by the 
interaction a A V. In the approximation of a slowly varying field 0, transitions in the 
system being studied can be considered as those that occur between the states \n) and 
which are induced by the interaction a 4 V. Since this interaction modifies the physical 
properties of the system, a finite number of its levels and their nonzero widths must be 
taken into account in calculating the probabilities W nm of the m(To) — > n(T) transitions. 
As a result, we arrive at 

W nm * |(n|W,(T,T )|m)| 2 exp{-r n AT}, (20) 

where AT = T — T and Ui is the evolution operator in the interaction representation 
[20]. Adiabaticity in the field enables us to consider specific transitions in the time 
interval AT that correspond to a given value of V(4>). 

The figure displays the total probability of universe decay, Wd ec = 1 — (Woo + Wo), 
and the quantity Wiq as calculated at V — 0.03, in which case there are only two levels 
in the system. It can be seen that, over the time interval AT < 50, the transitions in 
the system predominate and only for AT ~ 100 the probability that the universe tunnels 
through the barrier becomes commensurate with the probability that it undergoes the 
— > 1 transition in the prebarrier region. 

Since the rate at which the level width T n tends to zero is greater than the rate at 
which the potential decreases, the reduction of V with time results in that transitions 
become much more probable than tunnel decays, in which case the former fully determine 
the evolution of the quantum universe in the prebarrier region. If the universe has not 
tunneled through the barrier before the potential V of the field decreases to a value less 
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AT 



100 



Figure 1: Probabilities W\q and Wd ec versus the time interval AT = T — To at the 
parameter values of V = 0.03, ej, 0) = 3, ef } = 7, T = 2 x 10~ 6 , and r\ = 10~ 2 . 



than 0.01, a sufficiently large number of levels such that the probabilities of decays from 
them can be neglected are formed in it. Assuming that the amplitudes of transitions over 
the time interval AT are small, \V (n\a 4 \m}\ AT <C 1, we then find that 

Wn-i, n ' W n+2 ,n ' W n - 2 ,n ' W n+2 , n 

that is, the n — > n + 1 transition is more probable than the n — > n — 1 and n — > n + 2 
transitions. This means that the quantum universe can undergo transitions to ever 
higher levels with a nonzero probability. It is well known that the oscillator amplitude 
is quantized according to the condition a ~ \fn\ therefore, it can be concluded that the 
characteristic size a of the universe that did not undergo a tunnel transition increases as 
it is excited to higher levels. 

3.4- Parameters of the Early Universe 

In the adiabatic approximation, the expectation value a for the universe occurring in 
the lowest state y? is given by 



a w (y?oHy?o) = 



7T 



— V + O (V 2 ) 
16 v ; 



(22) 



whence it follows that a ~ 0.9 x 10 -33 cm for < V < 0.08. The value a determines the 
mean amplitude of oscillations of the classical universe filled with matter and radiation, 
thereby specifying its actual linear dimension. The maximal proper distance in a closed 
universe can be estimated atrfps7raps3x 10~ 33 cm; that is, the universe in the lowest 
state has a proper dimension on the order of the Planck length [19]. The presence of the 
minimal length removes the problem of the initial cosmological singularity. 

Equations (7) and (8), which are obtained within the general theory of relativity, 
also dictate the relationship between the quantities a, e and V. By using the value of 
e ~ 2.6, which we found for V ~ 0.08, we can estimate the classical turning points at 
di ~ 1.4 x 10 _33 cm and a 2 — 2.2 x 10~ 33 cm. The value of a x determines the maximal 
dimension of the universe occurring in the lowest state to the left of the barrier along 
the a axis, while a 2 characterizes its initial dimension upon tunneling from this state. In 
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this era, the matter- and radiation- energy densities are T® « V ~ 1.3 X 10 77 GeV/fm 3 , 
and T ° « 4& ~ 1.7 x 10 78 GeV/fm 3 , respectively; that is, we can see that, according 
to our model, the energy density in the early universe is determined primarily by the 
radiation-energy density. This result is fully consistent with what is commonly thought 
about the properties of the universe for a — > [9]. The total matter-energy density is 
p « 0.64 7«p, where m p is the Planck mass. Thus, we can see that quite reasonable results 
are obtained when the parameters e and V as derived on the basis of quantum theory 
are used in the equations of the general theory of relativity. 

It is interesting to estimate the quantity n at the a value coincident with the dimension 
of the presently observed part of the Universe. From the relation a ~ y/n at a ~ 10 28 cm, 
we obtain n ~ 10 122 , whence we can see that, if the quantum model being considered 
is extrapolated to the observed Universe, it occurs in a highly excited state. Quantum 
corrections to the classical equations of the general theory of relativity are extremely small 
in this case (they are of order ~ t/n). The resulting value of n ~ 10 122 is consistent with 
the estimates presented in [11, 21] and is confirmed by rigorous quantum-mechanical 
calculations within the radiation-free model that was considered in [19] and which is 
justified for the present, large, values of a at the matter-dominated stage. 



4. CONCLUSION 



The presence of radiation in the universe makes it possible to associate a privileged 
reference frame with it and to remove thereby an ambiguity in choosing the time coor- 
dinate. This opens new possibilities both in classical and in quantum cosmology. Upon 
performing quantization, there naturally arises the Schrodinger equation (13) with an 
effective interaction U in the form of a potential barrier. The evolution of the universe 
involves a quantum stage that is realized in the prebarrier region and which precedes the 
process of tunneling through the barrier. The dynamics of this stage is governed by the 
interaction of the gravitational and scalar fields. That the system in question possesses 
a spectrum of quantum (quasistationary) states and that transitions can occur between 
these states enable us to take a fresh look at the problem of the dimension of the Uni- 
verse. In the approach developed here, the universe is characterized by a minimal length, 
so that the singularity problem does not arise in it. The probability for the universe to 
undergo a tunnel transition is maximal in the lowest quantum state, where the energy 
density and the scale factor are on the same orders of magnitude as the corresponding 
Planck values. If a quantum universe tunnels from higher states, the dimensions of the 
region from which tunneling occurs can considerably exceed the Planck length. The con- 
stants E and V appearing in the Einstein equations are determined by the preceding, 
quantum stage. The use of the parameters in the general theory of relativity that were 
obtained on the basis of quantum theory leads to conclusions that are consistent with 
the currently adopted concepts of the early Universe and its subsequent evolution. 
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